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ABSTRACT 

It is well known that the accretion disk around a magnetized compact star can penetrate 
inside the magnetospheric boundary, so the magnetospheric radius Rq does not represent the 
true inner edge i?; n of the disk; but controversies exist in the literature concerning the relation 
between i?o and i?; n . In the model of Ghosh & Lamb, the width of the boundary layer is given 
by 5 = Rq — Rin Ciio, or R[ n ~ i?o, while Li & Wickramasinghe recently argued that R[ n could 
be significantly smaller than i?o in the case of a slow rotator. Here we show that if the star 
is able to absorb the angular momentum of disk plasma at Rq, appropriate for binary X-ray 
pulsars, the inner disk radius can be constrained by 0.8 ^ R m j Rq ^ 1, and the star reaches spin 
equilibrium with a relatively large value of the fastness parameter (~ 0.7 — 0.95). For accreting 
neutron stars in low-mass X-ray binaries (LMXBs), Rq is generally close to the stellar radius i?„ 
so that the toroidal field cannot transfer the spin-up torque efficiently to the star. In this case 
the critical fastness parameter becomes smaller, but R la is still near R . 

Subject headings: accretion, accretion disks - stars: magnetic fields -stars: neutron - stars: 
rotation 



1. Introduction 

Interaction between a magnetized, rotating star and a surrounding accretion disk has been of continued 
interest since the discovery of X-ray pulsars in binary systems (Giacconi et al. 1971, etc). Its solution 
may hold the key to understanding fundamental problems such as structure and dynamics of magnetized 
accretion disks, angular momentum exchange between the disk and the star, etc, in a variety of astrophysical 
systems, including cataclysmic variables, X-ray binaries, T Tauri stars and active galactic nuclei. A detail 
model for the disk-star interaction was first developed by Ghosh & Lamb (1979a, 1979b), who argued 
that the stellar magnetic field penetrates the disk via Kelvin-Helmholtz instability, turbulent diffusion and 
reconnection, producing a broad transition zone joining the unperturbed disk flow far from the star to the 
magnetospheric flow near the star. The transition zone is composed of two qualitatively different regions, a 
broad outer zone where the angular velocity is nearly Keplerian and a narrow inner zone or boundary layer 
where it departs significantly from the Keplerian value. The total torque N exerted on the star with mass 
M is then divided into two separate components: 

N = N in + N out , (1) 

where the subscript represents the appropriate zone in the disk. The "material" torque of the (inner) 
accretion flow is 

N in ~ A = MR 2 n k0 , (2) 
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where Slko = (GM / Rq 2 ) is the angular velocity of disk plasma at the magnetospheric radius Rq (G is the 
gravitation constant), and M is the accretion rate. Here we have adopted cylindrical coordinate (i?, <J>, z) 



centered on the star and the disk is assumed to be located on the z = plane, which is perpendicular to 
the star's spin and magnetic axes. The torque in the outer zone supplied by the Maxwell stress is given by 



where is the azimuthal field component generated by shear motion between the disk and the vertical 
field component B z . The torque N aut can be positive and negative depending on the value of the "fastness 
parameter" ui B = f2 s /f2ko = (Rq/R c ) 3 ^ 2 , where Q s is the angular velocity of the star and R c = (GM /Vt 2 ) 1 ^ 
is the corotation radius. 

Inside the magnetospheric boundary, the magnetic stress dominates the rotation of disk plasma, and 
forces it to corotate with the magnetosphere at the radius -RinQ To understand how far the disk can extend 
into the stellar magnetosphere and how the torque is transferred to the star, it is necessary to obtain 
the solutions of the structure of the magnetosphere. Ghosh & Lamb (1979a) suggested that Rq and S be 
constrained by the following relation 



(the subscript denotes quantities evaluated at R = Rq), which follows the equation of angular momentum 
conservation, though contribution from gas funnel flow to angular momentum transfer has not been 
included. Ghosh & Lamb (1979a) also assumed a priori that S corresponds to the electromagnetic screening 
length, and is much smaller than Rq. To testify the self-consistency of this assumption, from equation 
(4) we derive S ~ (Rq/ Ra) 7 ^ 2 Ro when lo s <C 1, where Ra is the conventional Alfven radius in spherical 
accretion. If Rq ~ 0.5i?A as obtained by Ghosh & Lamb (1979a), then indeed 5 <C Ro- However, if Rq ~ Ra 
as suggested by some recent investigations (e.g., Arons 1993; Ostriker & Shu 1995; Wang 1996; Li 1997), 5 
will be as large as Ro. 

Li, Wickramasinghe & Riidiger (1996) have presented a class of funnel flow solutions in which the 
angular momentum is carried by the matter. These authors argue that for an ideally conducting star, the 
stellar boundary condition requires the toroidal field be very small and most of the angular momentum of 
the matter in the funnel be propagated back to the disk rather to the star. Li & Wickramasinghe (1997) 
show that a consequence of this result is that the true inner disk radius i?i n could be much smaller than 
Rq, with R- m / Rq as low as ~ 0.1 — 0.2 in slow rotator case, and the values of the critical fastness parameter 
at which the net torque vanishes would also change. This "torqueless" suggestion was criticized by Wang 
(1997), who points out that a small toroidal field is enough to produce a significant torque on the star, and 
that the viscous stress is unable to remove angular momentum back outwards from the inner edge of the 
accretion disk. 

This paper is organized as follows. In section 2, based on Wang's (1997) arguments, we calculate the 
lower limit of R m that is required for efficient angular momentum transfer, and find it to be at least ~ 0.8i?o, 
implying that Rq is a good indicator of the inner edge of the disk. In section 3 we discuss the possibility 
of inefficient angular momentum transfer, which is likely to occur in LMXBs, and the implications for the 
resulting variation in the fastness parameter. We conclude in section 4. 




(3) 



{RB z B^/4n)o2TrRo26 ~ M(i?Au ) o 



(4) 



1 According to Spruit & Taam (1990), the accreting gas, nearly corotating with the star, could drift further inward across 
the field by an interchange instability. 
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2. The inner disk radius 

An estimation of i? can be obtained by setting the magnetic stress equal to the rate at which angular 
momentum is removed in the disk (Wang 1987) 

(i? 2 B,B z )o = [M^^]o, (5) 

where il = fl(R) is the angular velocity of disk plasma. Following Li & Wang (1996), we assume that Q(i?) 
reaches its maximum at R = Rq, i.e, (d£l/di?)o = 0, with fl(Ro) — f^ko, and rewrite equation (5) as 

R^qR-lq 2 
M(GMi? ) 1/2 R 3 ' 

The torque N- ln in equation (1) in fact contains two components: the real material torque Nf in the 
funnel flow, and the magnetospheric torque N mag arising from the shearing motion between the corotating 
magnetosphere and the non-Keplerian disk boundary layer, i.e., 

N in = Nf + A mag 

f-Ro ,-, t-Ro 

^R 2 ndR- B^B^dR, (7) 

where k — R m /Ro- 



kR$ ^ kR@ 



We first consider the second term on the righ-hand side of equation (7). In estimating the toroidal 
magnetic filed B^, we assume as usual that B^ is generated by the shearing between the disk and the stellar 
magnetosphere, and the growth of B^ is limited by diffusive decay due to turbulent mixing within the disk 
(e.g., Campbell 1992; Yi 1995; Wang 1995) 

b^ n s - n 

Several different phcnomcnological descriptions of disk-magnctosphere interaction (Wang 1995) cause little 
practical differences to our conclusions. Assuming that B z follows a dipolar field and using equations (6) 
and (8), we may write the magnetospheric spin- up torque into dimcnsionless form as 

^ mag = %^ = T s ii-(n s /n)} y -*dy, (9) 

Mo 0(1 - Ul s ) Jfe3/2 Ws 

where y = {R/R c f/ 2 . 

In the first term on the righ-hand side of equation (7), denotes the vertical mass loss rate from the 
boundary layer. This can be described by the equation (e.g., Ghosh & Lamb 1979a) 

Ifjf = 4nRpc s g(R), (10) 

scaling the vertical flow velocity from the boundary layer in terms of the local sound speed c s and 
introducing a "gate" function g(R) describing the radial profile of mass loss out of the boundary layer 
(p is the mass density). To keep things simple, we assume g(R) = outside Ro, and g(R) = 1 when 
Rin < R < Ro- Assuming that the disk plasma around Ro is thermally supported and optically thick to 
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Thomson scattering (appropriate for a binary X-ray pulsar) , from the a-disk model of Shakura & Sunyaev 
(1973), we obtain ^jf oc R' 1 or ^M. = { _^ k)R , which leads to 

Nf i n r 

^ = w = r^J kRo ^R7 dR - 

For generality, we have also assumed — const, and found that the results do not change significantly. 

The value of k can be obtained by setting £ = £f + £mag = 1. From equations (9) and (11) it can be 
seen that the magnitude of k depends on the detail profile of Ct(R) between R ln and Rq, which is unknown 
before solving the magnetohydrodynamics (MHD) equations. A lower limit of k, however, can be estimated 
in the following way. Because of the steep radial dependence of the stellar magnetic field (B z oc i?~ 3 ), the 
slope of the angular velocity ^ should increase from at Rq to a large value at R[ n . This requires that, in 
the Q(R)vs.R diagram (Figure la) Q(R) should always lie beyond the dashed line which represents a linear 
relation between Q, and R. As we know from equations (9) and (11), the larger f2, the larger the £, and the 
smaller the S required. So the smallest k corresponds to the extreme case that Q(R) decreases linearly from 
flko to f2 s between R and i?i n , i.e., 

n(R)=aR + b, (12) 

where 

O k0 (1 - UJ S ) 

a = lhJi~W' (13) 

and 

<»» 

In Figure lb we plot the calculated values of A: as a function of the fastness parameter w s in the dashed 
curve. It is seen that k naturally reaches 1 when lu s = 1, but remains ~ 0.8 within a large range of u) s . The 
solid curves in Figure la and lb represent an example with more practical boundary conditions, which are 
always located above the dashed curves. 



3. The fastness parameter 

Our calculations in last section are based on the assumption that a torque of order Mi?o 2 ^ko can be 
efficiently transmited to the star by both the gas funnel flow and the magnetic stress. It has been suggested 
that (Shu et al. 1994; Ostriker & Shu 1995; Li, Wickramasinghe, & Riidiger 1996) a star may accrete 
matter from a magnetically truncated Keplerian disk without experiencing any spin-up torque, because it 
would require the surface magnetic field to have a large azimuthal twist which would lead to dynamical 
instabilities. Wang (1997) estimated the azimuthal pitch 7* at the stellar surface required to transmit this 
torque to be 

^^J^R^s/zRo r ._x 
7 *~ B n ~ 23/2 ( i? j 5' (i5j 

where the subscript asterisk denotes quantities at the stellar surface. Since MHD instabilities are expected 
to occur only when |7*| £ 1, one sees that, with S ^ 0.2i? and £ ~ 1, 7, could be much less than unity if 
i?* <C Rq- This condition is satisfied by binary X-ray pulsars, which generally possess surface magnetic fields 
as strong as 10 12 — 10 13 G, and hence Rq ~ 10 8 — 10 9 cm with the mass accretion rates of 10 16 — 10 18 gs _1 , 
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Fig. 1. — (a) The radial profile of the angular velocity £l(r) of disk plasma in the boundary layer. The 
dashed line denotes a linear decrease of Q(R) with R, Q — R relations above which are physically allowed. 
The solid curve represents one of the practical cases, (b) The radius ratio of k = R m /R$ as a function of 
the fastness parameter u B . The solid and dashed curves correspond to those in (a), respectively. The latter 
puts a lower limit on k. 

much larger than i?„ ~ 10 6 cm. Thus the standard accretion disk models based on Ghosh & Lamb (1979a, 
1979b) are suitable for these systems, and the value of the critical fastness parameter uj c , according to Wang 
(1995) and Li & Wang (1996), may lie between - 0.7 and 0.95. 

There could exist possibility of an accreting star unable to absorb the excess angular momentum at 
i?o- This may occur in accreting neutron stars in LMXBs, in which the condition |7*| -C 1 may not hold 
unless £ ~ 0, because the surface field strengths in these neutron stars are so low (10 8 — 10 9 G or lower) 
that the accretion disk usually extends close to the star's surface, i.e., R ~ i?*. Since MHD instabilies 
prevent transmission of the torque of -/V in , in this case the star only experiences the magnetic torque N out , 
and the critical fastness parameter decreases to ~ 0.4 — 0.6 (Li & Wickramasinghc 1997). Note that the 
angular momentum of the material at Rq cannot be transported from the boundary layer to the outer parts 
of the disk by either viscous stress or a stellar magnetic field threading the inner region of the disk (Wang 
1997), it is more likely to be removed by a magnetocentrifugal wind. During wind mass loss, some angular 
momentum of the star could also be lost, equivalent to an extra braking torque exerted on the star, so that 
the critical fastness parameter would be even smaller. 

Observational evidence for the above arguments comes from kiloHertz quasi-periodic oscillations (kHz 
QPOs) recently discovered with the Rossi X-ray Timing Explorer (RXTE) in some 16 neutron star LMXBs 
(see van der Klis 1997 for a review). These kHz QPOs are characterized by their high levels of coherence 
(with quality factors Q = v/Av up to ~ 200), large rms amplitudes (up to several 10%), and wide span of 
frequencies (~ 300 — 1200 Hz), which, in most cases, are strongly correlated with X-ray fluxes. In many 
cases, two simultaneous kHz peaks are observed in the power spectra of the X-ray count rate variations, 
with the separation frequency roughly constant (e.g., Strohmayer et al. 1996; Ford et al. 1997; Wijnands 
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et al. 1997). Sometimes a third kHz peak is detected in a few atoll sources during type I X-ray bursts at 
a frequency equal to the separation frequency of the two peaks (Strohmayer et al. 1996; Ford et al. 1997) 
or twice that (Wijnands & van der Klis 1997; Smith, Morgan & Bradt 1997). This strongly suggests a 
beat-frequency interpretation, with the third peak at the neutron star spin frequency (or twice that), the 
upper kHz peak at the Keplerian orbital frequency at a preferred radius around the neutron star, and the 
lower kHz peak at the difference frequency between themQ. If the upper kHz QPOs are associated with 
the orbital frequency at the magnetospheric radius or the sonic point in the accretion disk (Strohmayer et 
al. 1996; Miller et al. 1998; see also Lai 1998), one would expect that the disk can extend towards the 
marginally stable orbit of a neutron star. 

As pointed by Zhang et al. (1997) and White & Zhang (1997), these neutron stars may have accreted 
a considerable amount of mass, so they possibly have reached the equilibrium spin. The measured spin and 
QPO frequencies of these neutron stars reveal the critical fastness parameter uj c lying between ~ 0.2 and 
0.7. Additionally, considering the spin and orbital evolution in low-mass binary pulsars which are thought 
to originate from LMXBs, Burderi et al. (1996) also found a small value of uo c ~ 0.1 from the observed 
relation between spin period, magnetic field and orbital period. Though subject to both observational and 
theoretical uncertainties, these results do indicate a smaller fastness parameter than predicted by current 
theory for binary X-ray pulsars. 

4. Conclusions 

We summarize our results as follows: (1) The inner disk radius R\ n is always close to the magnetospheric 
radius Rq. (2) For binary X-ray pulsars, the critical fastness parameter lies in the rang 0.71 — 0.95; neutron 
stars in LMXBs may not absorb angular momentum flux at Rq as efficiently as X-ray pulsars, resulting in a 
critical fastness parameter considerably smaller than in X-ray pulsars. 

This work was supported by National Natural Science Foundation of China. 



2 The applicability of such a model is less clear when the difference frequency is not constant as in the Z source Sco X-l and 
possibly the atoll source 4U 1608-52; see van der Klis et al. (1996) and Mendez et al. (1998). 
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